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Here, we present a detailed study on the temperature and angular dependence of the Shubnikov- 
de-Haas (SdH) effect in the semi-metal WTe 2 . This compound was recently shown to display a 
very large non-saturating magnetoresistance which was attributed to nearly perfectly compensated 
densities of electrons and holes. We observe four fundamental SdH frequencies and attribute them 
to spin-orbit split, electron- and hole-like, Fermi surface (FS) cross-sectional areas. Their angular 
dependence is consistent with ellipsoidal FSs with volumes implying a modest excess in the density 
of electrons with respect to that of the holes. We show that density functional theory (DFT) calcu¬ 
lations fail to correctly describe the FSs of WTe 2 . When their cross-sectional areas are adjusted to 
reflect the experimental data, the resulting volumes of the electron/hole FSs obtained from the DFT 
calculations would imply a pronounced imbalance between the densities of electrons and holes. We 
find evidence for field-dependent Fermi surface cross-sectional areas by fitting the oscillatory com¬ 
ponent superimposed onto the magnetoresistivity signal to several Lifshitz-Kosevich components. 

We also observe a pronounced field-induced renormalization of the effective masses. Taken together, 
our observations suggest that the electronic structure of WTe 2 evolves with the magnetic field. This 
evolution might be a factor contributing to its pronounced magnetoresistivity. 


I. INTRODUCTION 

Transition metal dichalcogenides TMX 2 (where TM 
stands for a transition metal, and X for a chalcogen) have 
become the subject of intense research activity since, as 
graphene, they can be exfoliated down to a single atomic 
layer [l|. Semiconducting compounds present remark¬ 
able optoelectronic properties [2] while the whole family 
displays a wide range of electronic states [1]; from wide 
gap semiconductors such as trigonal 1 T-HfS 2 3], to band 
gaps comparable to Si such as in <a-MoTe 2 [J], or Mott 
systems such as lT-TaS 2 [5j]. Some compounds display 
charge-density wave(s) that coexist with superconductiv¬ 
ity as in iF-NbSe 2 @, while others are semi-metals such 
as /3-MoTe 2 or WTe 2 0. 

WTe 2 , which crystallizes in a distorted, orthorhom¬ 
bic variant of the Cdl 2 structure with octahedral co¬ 
ordination around W, was recently reported to display 
an extremely pronounced and non-saturating magnetore¬ 
sistance which [7]. Following simple arguments already 
discussed by PippardQ, this effect was attributed to a 
nearly perfect compensation between the density of elec¬ 
trons and holes. Evidence for electron/hole compensa¬ 
tion was recently provided by angle-resolved photoemis¬ 
sion (ARPES) experiments claiming to observe electron- 
and hole-like Fermi surface pockets of almost the same 
size [ 9 ]. These results contrast markedly with a subse¬ 
quent ARPES study [l0] which claim to observe up to 
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9 FS sheets, including a hole-pocket at the Brillouin 
zone center or the T-point, and two hole-pockets and two 
electron-pockets on each side of T and along the T-X di¬ 
rection. The same study finds evidence for spin-orbit 
split bands and for an opposite texture between the spin 
and the angular orbital momentum around the T-point 
which should suppress carrier backscatteringfioj. The 
application of a magnetic field would not only alter the 
spin texture, it would also displace, through the Zeeman 
effect, the spin-split bands relative to the Fermi level, 
opening backscattering channels and hence increasing the 
resistivity. Similar arguments were invoked to explain the 
pronounced magnetoresistivity observed in CdsAs^Tl]. 

Electrical transport studies focused on the Shubnikov- 
de Haas-effect yielding results UMi which contrast 
markedly with both ARPES studies. The first study re¬ 
ports the observation of just 3 Fermi surface pockets in 
marked contrast with band structure calculations which 
alwayspredict an even number of FS pockets for bulk 
WTe 2 [9|, [ 13 , [fij. Subsequent studies, more in line with 
the calculations, observe 4 Fermi surface pockets, or two 
sets of concentric electron- and hole-like FSs[l3l-[l5| at ei¬ 
ther side of the T— point. The angular dependence of the 
observed frequencies, or of the FS cross-sectional areas 
obtained from the Onsager relation, suggests ellipsoidal 
pockets with volumes yielding an excess of 4 % in the 
density of holes when compared that of the electrons. 
This imbalance would be comparable to that of elemen¬ 
tal Bi which, in contrast to WTe 2 , displays a saturating 
magnetoresistance jl7] . 

Subsequent calculations H8| confirm the previously re¬ 
ported FS geometry @,0111, while analyzing the topo- 
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logical properties of its electronic structure around the 
Fermi level. Ref. finds a linear touching between hole 
and electron Fermi-surfaces leading to Berry phase cur¬ 
vature singularities, or to two sets of four Wey 1-points in 
the first Brillouin zone, each Wey 1-point characterized by 
topological charges of ±1. At the moment it remains un¬ 
clear how such electronic structure would contribute to 
WTe 2 pronounced magnetoresistivity. Ref. m also finds 
that its FS is particularly sensitive to the exact position 
of the Fermi level with small displacements in sp 
leading to topological Lifshitz transitions, in agreement 
with our DFT calculations being discussed within this 
manuscript. 

Here, we present angle- and temperature-dependent 
magnetoresistivity measurements in fields up to H = 35 
T on WTe 2 single-crystals, in an attempt to elucidate the 
number of Fermi surface sheets. In parallel, we performed 
density-functional theory (DFT) calculations incorporat¬ 
ing the spin-orbit effect. By adjusting the position of 
the chemical potential we attempted to precisely describe 
the geometry of the measured Fermi surfaces. In agree¬ 
ment with a previous report[l3] we observe four funda¬ 
mental Shubnikov-de-Haas frequencies, or Fermi surface 
cross-sectional areas, which we ascribe to spin-orbit split 
electron- and hole-pockets. As argued in Ref. 0, their 
angular dependence can be described in terms of ellip¬ 
soidal FSs. Therefore, we find that the DFT calculations 
are unable to correctly describe the angular dependence 
of the FS cross-sectional areas, particularly in what con¬ 
cerns the hole-pockets. For fields aligned along the c— 
and the b —axis we observe the effective masses to be¬ 
come renormalized, i.e. considerably heavier as the field 
increases. As the FFT spectrum of the oscillatory signal, 
taken on a limited field window, is unable to resolve all 
cross-sectional areas, we fit it to five and four Lifshitz- 
Kosevich components. Our fits clearly indicate that one 
cannot fit the SdH signal over the entire field range by 
maintaining the frequencies and the effective masses as 
constants. Good fits are obtained within a limited field 
range, confirming that the effective masses (and likely 
the cross-sectional areas) are field-dependent. 


II. EXPERIMENTAL 

WTe 2 was grown via a chemical vapor transport tech¬ 
nique using chlorine as the transport agent: W (99.999 
%) and Te (99.9999 %) were heated up to a peak temper¬ 
ature of 750 °C at a rate of 100 °C/h, and subsequently 
cooled down to room temperature at a rate of 10 °C/h. 
The obtained powder was ground with a mortar and pes¬ 
tle and heated again following the the same temperature 
profile. The remaining product was subsequently com¬ 
bined with TeCU (99.999 %, 1.5 mg/cm 3 ) and heated up 
to a peak temperature of 750 °C at a rate of 100 °C/h 
under a temperature gradient of ~ 50 °C. The crystals 
used for this study displayed resistivity ratios in excess of 
10 2 . Subsequent crystals grown via the Te flux method, 


yielded much higher RRRs (~ 700) and show similar be¬ 
havior (see supplementary information). 

The synthesis procedure yielded platelet like single 
crystals, several millimeters in length and typically a mil¬ 
limeter wide, making the a-, 6-, and the c-axis easily 
distinguishable. The stoichiometric composition was de¬ 
termined by energy dispersive X-ray spectroscopy and 
single-crystal X-ray refinement. A conventional four ter¬ 
minal configuration was used for the resistivity measure¬ 
ments which were performed under high magnetic-fields 
either by using a 35 T resistive magnet or a Physical 
Property Measurement System (PPMS). 
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FIG. 1. (a)Resistivity p as a function of the magnetic field 
H applied along the c-axis, for a WTe 2 single-crystal and for 
three representative values of the temperature T, i.e. 0.5 K 
(black line), 6 K (blue) and 10 K (red), respectively. Here, 
the current is applied along the a-axis, hence perpendicularly 
to H. (b)Same as in (a) but for fields along the 6-axis of 
the crystal, (c) Same as in (a) but for currents parallel to 
H which is applied along the a-axis. (d) Shubnikov-de-Haas 
signal superimposed onto the curves in (a) after subtraction 
of a polynomial background, (e) Same as in (d) but from the 
curves in (b). (f) Same as in (d) but from the curves in (c). 

Figures [1] (a), (b), and (c) show the resistivity p of 
a WTe 2 single-crystal as a function of the magnetic 
field applied along the a, 6, and c-axis respectively, 
for several temperatures and with the electrical cur¬ 
rent applied along the a-axis. In agreement with pre- 
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vious reports[7, [HI, fields along the c-axis yield the 
largest magnetoresistivity, although in our crystals the 
anisotropy in magnetoresistivity between fields along the 
c- and the 6-axis is jus t a factor of 4 and not beyond one 
order of magnitude[13|- The smallest change in resistiv¬ 
ity is observed for fields applied along the a-axis, since in 
this configuration the Lorentz force is expected to vanish. 
For all three sets of data the oscillatory component super¬ 
imposed onto the smoothly increasing background corre¬ 
sponds to the Shubnikov-de-Haas effect. Figures [1] (d), 
(e), (f), display the oscillatory component superimposed 
onto p(H,T) as a function of the inverse field iJ _1 , for 
fields oriented along all three crystallographic axes and 
for several temperatures (from the curves in (a), (b) and 
(c), respectively). To obtain the oscillatory component, 
the smoothly varying background was fit to a polynomial 
which was subsequently subtracted from the raw data. 

Figures [2] (a), (b), and (c) show the fast Fourier trans¬ 
form (FFT) of the SdH oscillatory signal displayed in 
Figs. Q](d), (e), and (f), respectively. As seen, for fields 
along all three crystallographic axes four fundamental 
frequencies are clearly observed. This is in agreement 
with the observations of Ref. [l3| but contrasts markedly 
with both ARPES studies H and JlO. For fields along 
the c-axis one observes peak frequencies at F a i = 90 T, 
Fpi = 130 T, Fp 2 = 145 T, and F a2 = 160 T, respec¬ 
tively. As we discuss below, and based on DFT calcu¬ 
lations, the a and the /3 frequencies are identified with 
spin-orbit split hole- and electron-pockets, respectively. 
The indexes 1 and 2 distinguish the smaller from the 
larger spin-orbit split Fermi surfaces. The DFT calcu¬ 
lations indicate the presence of identical sets of electron 
and hole pockets at either side of the T—point making 
it a total of 8 Fermi surface sheets in the first Brillouin 
zone. We see a gradual increase in the frequencies as the 
crystal is rotated from c- towards the 6-axis or the a-axis. 
For fields along the 6-axis the peaks are shifted towards 
F a i = 180 T, Fpi = 220 T, Fp 2 = 236 T, and F a2 = 262 
T, respectively. For fields along the a-axis the peaks are 
observed at F a i = 250 T, Fpi = 300 T, Fp 2 = 357 T, 
and F a2 = 390 T, respectively. Ref. claims the exis¬ 
tence of 9 Fermi surface pockets in contrast to the eight 
pockets implied by our experiments when combined with 
the calculations described below. Here, we just mention 
that another very small electron pocket do emerge at the 
T—point by shifting the Fermi level in the DFT calcula¬ 
tions with respect to the position required to match our 
observations. Indeed this would lead to a total number of 
9 Fermi surface sheets. Nevertheless, their cross-sectional 
areas would strongly disagree with our experimental ob¬ 
servations. Furthermore, the different ARPES studies 
SGI disagree among themselves on the existence of this 
extra pocket. As we show below, the DFT calculations 
indicate that the geometry of the Fermi surface(s) is ex¬ 
tremely sensitive with respect to small displacements of 
the Fermi level, i.e. of the order of 20 meV, which is com¬ 
parable to the experimental resolution of ARPES. The 
distinction among the ARPES results is attributable to 
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FIG. 2. (a) Fast Fourier transform of the SdH signal super¬ 
imposed onto p(H) for fields along the a-axis and for T — 1.3 
K. Here, the oscillatory signal was extracted from a field in¬ 
terval ranging from 3 to 35 T. As is discussed in the main 
text, the indexation of the four peaks observed in the FFT 
spectrum is based upon band structure calculations, where F a 
and Fp denotes hole- and electron-pockets respectively, and 
the indexes 1 and 2 differentiate spin-orbit split Fermi surface 
pockets, (b) Same as in (a) but for fields along the 6-axis and 
for T = 0.5 K. (c) Same as in (a) but for fields along the 
c-axis and for T — 0.5 K. Inset: FFT spectra extracted for 
different field intervals. Notice the loss of frequency resolu¬ 
tion when the FFT is extracted either from traces acquired 
under high fields or from a limited field range. To allow the 
visual comparison between the different FFT spectra the blue 
and black traces were multiplied by a factor 75 and 3.5, re¬ 
spectively. Notice how within the interval AH = 8 — 20 T 
the peak corresponding to the magnetic breakdown orbit F$ 
becomes more pronounced than the fundamental peaks. 


the resolution between the different set-ups. 

For fields along the c-axis one also observes a pro¬ 
nounced fifth peak, F$ ~ 250 T, which clearly corre¬ 
sponds to F a i + F a 2 . This frequency was previously 
attributed to magnetic breakdown or to the magnetic 
field-induced carrier tunneling between both orbits[l3j- 
However, the fact that the amplitude of the peak associ¬ 
ated with F§ is more pronounced than peaks associated 
with the fundamental orbits, such as F a2 , contradicts this 
scenario, unless both FS sheets were in close proximity 
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FIG. 3. (a) Magnitude of the FFT peaks extracted from the oscillatory component superimposed into the magnetoresistivity 
measured for fields applied along the a-axis and as a function of the temperature T. Solid lines are linear fits to the Lifshitz- 
Kosevich formalism from which we extract the effective masses /i. (b) Same as in (a) but for fields along the 6-axis, (c) Same as 
in (a) but for fields along the c-axis. (d) Magnitude of FFT peak observed at F§ which corresponds to a magnetic breakdown 
orbit. 


at the point of nearly touching. In fact, as seen in the 
inset of Fig. [2] (c), which displays the FFT spectra ac¬ 
quired under distinct field intervals, this orbit becomes 
progressively more pronounced at higher fields, which is 
consistent with magnetic breakdown. However, as illus¬ 
trated by the FFT spectrum extracted from the field in¬ 
terval ranging from 8 to 20 T, the amplitude of F$ can, in 
some field intervals, become more pronounced than the 
amplitude of all other fundamental peaks. This suggests 
that the Zeeman-effect is altering the geometry of the 
<a-pockets, hence decreasing their (spin-orbit induced) 
splitting in some region(s) of &-space. Our DFT anal¬ 
ysis, discussed below, is indeed consistent with this pos¬ 
sibility. This would imply that the electronic structure 
at the Fermi level evolves with field due to the Zeeman- 
effect. We considered this hypothesis and evaluated the 
possibility of studying in detail the evolution of the SdH 
frequencies, and of the concomitant effective masses, as 
a function of the field. But as illustrated by the inset of 
Fig. [2] (c), the resolution among peaks in the FFT spec¬ 
tra, or the peak sharpness, decreases considerably when 
the very low field interval is excluded from the oscillatory 
signal, thus preventing such a detailed analysis. 

Figures [3] (a), (b), and (c) display the amplitude of 
each peak observed in the FFT spectra as a function of 
the temperature T, for fields along the a—, 6—, and the 
c— axis, respectively. Solid lines are fits to the Lifshitz- 


Kosevich formula A/ sinhA (where X = 27t 2 &bT /huj c = 
14.69 fiT/H, with /i being the effective mass (i and H the 
average field value), from which we extract the effective 
masses. As seen, the fits are excellent and yield pre¬ 
cise values for (i. More importantly, according to the fits 
in Fig. [3] (d), the effective mass of or //<$, is found to 
evolve as a function of the field: i.e. /is = (0.73±0.05)rao 
for F[ ranging from 3 to 13 T, fi# = (0.9 db 0.07)rao for 
fields between 8 and 20 T, ns = (0.97d=0.02)rao for fields 
ranging from 15 to 35 T. In Table I below summarizes 
the experimental information extracted from this study, 
including effective masses for H || c— and to the b— axis, 
in each case extracted from two intervals in field, i.e. a 
low field region and the whole field window up to 35 T. 
As seen, the effective masses extracted from the whole 
field range are approximately a factor of 2 heavier than 
those extracted from the low field region, clearly indicat¬ 
ing a renormalization of the effective mass induced by 
the application of a magnetic field. In the same table 
we include the results of an attempt to match the experi¬ 
mentally determined FSs with those determined from the 
DFT calculations by independently adjusting the posi¬ 
tion of the electron- and hole-bands relative to the Fermi 
level. 

A renormalization of the effective masses necessarily 
implies the evolution of the electronic structure at the 
Fermi level due to the application of an external field. 
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FIG. 4. (a) Kohler plot of the resistivity p, normalized by its zero field value po, as a function of (H / po ) 2 for fields applied 
along the a—axis, (b) Same as in (a) but for fields along the 6-axis, (c) Same as in (a) but for fields along the c-axis. 


Such an evolution is likely to affect the geometry of the 
Fermi surface. To test this hypothesis in Fig. [TT] we 
attempt to fit the oscillatory component, or the SdH 
signal superimposed onto the magnetoresistivity to the 
five frequencies observed in the FFT spectrum. How¬ 
ever as illustrated by Figs. 5 (a) and (b) it is not pos¬ 
sible to fit the SdH signal over the entire field range 
to 5 Lifshitz-Kosevich components (4 fundamental fre¬ 
quencies plus Fs), each composed of a sinusoidal term 
whose amplitude is normalized by i7 -1 / 2 , the exponen¬ 
tial Dingle factor (whose argument is effective mass and 
scattering rate dependent) and the temperature factor 
oc [sinh (ap J T/H)]~ 1 , where a is a constant. As seen in 
Fig. 5 (a), for fields below H m 10 T one can fix the fre¬ 
quencies to the values obtained from the FFT spectrum 
and obtain a good fit of the experimental data over an 
extended range in H~ l if one considers the experimen¬ 
tal noise and a less than perfect background subtraction. 
However, as seen in Fig. 5 (b) if one attempts the same 
exercise for fields beyond 10 T one cannot converge to a 
fit describing the oscillatory signal over the entire H~ x 
range. This is particularly surprising given that this win¬ 
dow in inverse field AH -1 ~ 0.714 is shorter than the 
interval AH~ X > 0.16 previously used for Fig. 5(a). As 
illustrated by Figs. 5 (c) and (d) it is nevertheless possi¬ 
ble to converge to good fits if one splits this interval into 
two intervals, i.e. from 10 to 18.6 T and from 18.6 to 
35 T. For fields beyond 10 T we were unable to simulta¬ 
neously resolve F a 2 and F^\ the fits systematically con¬ 
verge to an average value between both frequencies and 
to another frequency around 400 T (curiously detected 
in the FFT spectra shown in the inset of Fig. 2(c)). 
This is likely attributable to the limited interval in H~ x . 
Hence, we have re-done the analysis after removing one 
LK component (blue traces) which does not compromise 
the fits. The important observation is that F a 1 increases 
from 93 T (in the field interval H ~ 3.5 to 10 T) to 
103 T (H ~ 18.6 to 35 T), F pi from 132 to 145 T and 
(F a2 + F> 2 )/2 = (150 + 164)/2 = 157 to 187 T. One could 
argue that the fits of the highest field oscillatory data are 


not excellent and hence question the validity of the ob¬ 
tained frequencies. However, here we have shown that 
it is not possible to fit the entire field range with fixed 
frequencies and effective masses: the Dingle and temper¬ 
ature factors also depend on the field interval chosen for 
the fits. Together with the previously shown effective 
mass renormalization, this indicates that the electronic 
structure of WTe 2 evolves with field due to the Zeeman 
splitting. 

In Fig. [4] we show Kohler plots of the resistivity p nor¬ 
malized by its zero field value po as a function of the 
square of the field H 2 normalized by po and for fields 
along all three crystallographic orientations. Recently, it 
was claimed[l4| that WTe 2 would not obey the Kohler 
rule [§J, leading to an anomalous linear magnetoresistance 
for currents parallel to the field. However, as seen in 
Fig. SJ our data clearly follows Kohler scaling with very 
mild deviations for fields along the b— axis, and only for 
the curves taken at higher temperatures. At lower fields, 
and for fields along the a—axis, indeed p(H) ~ H. Since 
for this current/field configuration the Lorentz force is 
nearly zero, which leads to no orbital magnetoresistiv¬ 
ity, a very careful experimental effort, such as very pre¬ 
cise field-alignment and homogeneous current distribu¬ 
tion will be required to characterize this effect. The fact 
that WTe 2 nearly obeys Kohler rule for all three crys¬ 
tallographic orientations implies that its orbital magne¬ 
toresistivity is dominated by scattering by phonons and 
impurities. 


III. DFT CALCULATIONS AND THE 
EVOLUTION OF THE FERMI SURFACES 

One of the main goals of this study is to compare the 
geometry of the experimentally detected FSs with those 
resulting from ab-initio calculations such as the ones in¬ 
cluded in Refs. 0 I 9 L fl3l and [l6|. Notice that the initial 
calculations 0, 0, [l6| ignored spin-orbit coupling, im¬ 
plying the existence of one hole- and one electron-like 
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FIG. 5. (a) Raw SdH signal (black trace) at T — 0.5 K as a function of the inverse field i/ _1 , for fields ranging between 
3.5 and 10 T applied along the c-axis. Magenta line corresponds to a fit of the SdH signal to 5 Lifshitz-Kosevich oscillatory 
components, each having one of the frequencies observed in the FFT spectrum, (b) Same as in (a) but for fields ranging from 
10 to 35 T. In contrast to what is observed in (a) the fit does not correctly describe the experimental data over the entire held 
range, (d) Same as in (b) but over a smaller held range, i.e. from 10 to 18.6 T. Blue line depicts a ht to a 4 LK components, 
where one of the components corresponds to an average value between F a 2 and Fp 2 - In contrast to what is seen in (b) both 
hts are able to reproduce the experimental data, (d) Same as in (c) but for helds ranging between 18.6 and 35 T. 


pocket. This would be at odds with the results reported 
by Ref. Qjl In this paper, electronic structure calcula¬ 
tions were performed by using the Vienna ab-initio sim¬ 
ulation package jl9l-[22j(VASP) within the generalized 
gradient approximation (GGA). We have included the 
contribution of spin-orbit coupling in our calculations. 
The Perdew-Burke-Ernzerhof (PBE) exchange correla¬ 
tion functional [23] and the projected augmented wave 
(PAW) methodology [24| was used to describe the core 
electrons. The 5 d and 6 s electrons for W and the 5s and 
5 p electrons for Te were treated as valence electrons in 
all our calculations. The energy cut off for the plane- 
wave basis was chosen to be 600 meV. A total of 96 
bands and a k —point mesh of 8 x 8 x 8 were used for 
the self-consistent ground state calculations. A total of 
100 points were chosen between each pair of special 
A:—points in the Brillouin-zone for the band-structure cal¬ 
culations. 

In order to obtain an accurate and detailed represen¬ 


tation of the Fermi surface we have used a 80 x 80 x 80 
k-point mesh. The areas of the observed and calculated 
FS orbits are related to the SdH frequencies by the On- 
sager relation, 


2tt 2 k 


( 1 ) 


where <po = 2.07 x 10 -15 Tm 2 is the quantum of flux and 
A £ is the area of the electron orbit perpendicular to the 
field. The effective masses of the electron- and the hole- 
pockets at these extremal orbits are determined by using 
the relation: 


m l (e) 


ft 2 dA s (e) 
2?r de 


( 2 ) 


We have calculated the areas of the orbits by consider¬ 
ing cuts of the Fermi surface with planes intersecting the 
Fermi surface perpendicularly to the a-, b- and the c-axis 
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FIG. 6. (a) Fermi surface of WTe 2 . The spin-orbit split elec¬ 
tron pockets are shown in blue, and in transparent light blue 
color, respectively. The hole pockets are shown in red and in 
pink color. The position of the energy of the electron and hole 
pockets was independently adjusted to obtain a good agree¬ 
ment with the experimentally observed Shubnikov de Haas 
frequencies for magnetic fields along the c-axis. In order to 
obtain high resolution, the size of the fc-point mesh used was 
80 x 80 x 80 . This is required to obtain, as accurately as possi¬ 
ble, the angle dependence of the frequencies for both electron 
and hole orbits. The reason for 2 electron- and 2 hole-like 
pockets, as opposed to just two doubly degenerate ones, is 
the spin-orbit coupling. Therefore, the spin-orbit coupling is 
particularly large (relative to the energy dispersion near the 
Fermi level) for the hole pockets, (b) The electron pockets 
are shown in a higher resolution, larger image, (c) Same as 
in (b) but for the hole-pockets. 


at various different angles. In an attempt to match the 
experimentally observed FSs with the calculated ones, 
we chose to adjust our calculations to the experimen¬ 
tal cross-sections obtained for fields along the c-axis. 
For this, the energy of the electron pockets were down 
shifted with respect to the Fermi level by an amount of 


25 meV, while the hole pockets were displaced upwards 
by 20 meV. The fermi surfaces were generated using the 
eigenvalues obtained from VASP and were visualized us¬ 
ing the XCrysden software pHj |. As seen in Fig. Mi the 
resulting Fermi surface is composed of two sets of spin- 
orbit split pockets at either side of the T-point and along 
the T-X direction, with the electron-like pockets being 
depicted in blue or clear blue color, and the hole-like in 
red or in clear and nearly transparent red. 

In Table [III below, we summarize the experimental re¬ 
sults as well as those resulting from the DFT calculations. 
In particular, we present a comparison between the ex¬ 
perimental effective masses /i, extracted from the fits in 
Fig. [3j and the band masses fi h resulting from the DFT 
calculations. Here, the important observation is not only 
the marked contrast between the band masses and the 
experimental effective-masses with the band masses be¬ 
ing systematically heavier (in particular for the <a-orbits), 
but also their field dependence. As seen in Table an the 
effective masses are strongly field dependent, increasing 
by nearly a factor of 2 when they are extracted from the 
whole field range as compared with the masses extracted 
from the low field region. This is a very clear indication 
for the role of the magnetic field in altering the electronic 
structure at the Fermi level and in leading to a concomi¬ 
tant increase in the density of states at the Fermi level 
q(Ef). For samples that start to display quantum oscil¬ 
latory behavior under fields H < 3 T it is remarkable, 
as seen in the inset of Fig. [2] (c), that one cannot resolve 
all four peaks in the FFT spectrum when it is taken over 
a broad range in i7, such asi7 = 10—)>35T. However, 
as illustrated in the Supplemental Material [26| through 
simple simulations of the experimental data, this just re¬ 
flects the much lower number of oscillations within this 
inverse field range when compared to the number of os¬ 
cillations observed in fields ranging from < 3 T up to 9 
T. 
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Me 
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35) T 

U/ellip. 

nDFT 

Foci 

1.37 

(1.17 

± 0.06) 

0.84 

(0.44 

± 

0.11) 

(0.97 

± 

0.04) 

0.72 

(0.42 

± 

0.01) 

(0.79 

± 

0.04) 

0.57 

0.78 

F a . 2 

1.32 

(1.20 

± 0.07) 

1.02 

(0.44 

± 

0.05) 

(1.13 

± 

0.05) 

1.1 

(0.47 

± 

0.02) 

(0.95 

± 

0.05) 

1.14 

2.31 

Fpi 

1.02 

(0.64 

±0.14) 

0.56 

(0.49 

lb 

0.01) 

(0.84 

± 

0.03) 

0.43 

(0.33 

± 

0.01) 

(0.62 

± 

0.03) 

0.84 

3.9 

Fp 2 

0.97 

(0.75 

± 0.04) 

0.57 

(0.58 

± 

0.03) 

(0.92 

± 

0.04) 

0.51 

(0.36 

± 

0.01) 

(0.64 

± 

0.04) 

1 

5.14 


TABLE I. Comparison between the experimental results and the results of the DFT-calculations. li h a ^, c are the band masses 
from the DFT calculations for orbits perpendicular to fields applied along all three crystallographic axes (in units of free electron 
mass mo), li a ,b,c are the experimental effective-masses, n e iii p . are the carrier densities as estimated from a simple ellipsoidal 
approximation of the Fermi surfaces, and tidft are the carrier densities resulting from the DFT calculations. The ns are 
expressed in units of 10 19 cm -3 . 


Table HU also shows the density of carriers n e iii P . that 
one would extract from our experiments by assuming el¬ 


lipsoidal FS pocketsjl3]. As discussed further below, we 
fit the observed angular dependence of the cross-sectional 
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FIG. 7. Position of the main four peaks detected in fast 
Fourier spectra as a function of the orientation of the field 
relative to the crystallographic axes. In order to compare 
experimental results with the theoretical calculations dashed 
lines depict the DFT calculated angular dependence of the 
frequencies for both electron and hole pockets. Notice the 
overall poor agreement with the calculations. Solid lines cor¬ 
respond to fits assuming ellipsoidal Fermi surfaces. 


areas to ellipsoidal Fermi surfaces to extract the corre¬ 
sponding Fermi-vectors k F b,c . The volume of each FS 
pocket is estimated through the expression for the vol¬ 
ume of an ellipsoid, i.e. Vps = 4:7rk F k F k F /3. To extract 
rieiiip. one renormalizes by the volume of the Brillouin 
zone. This simple estimation yields a ^ 7% larger density 
of electrons with respect to that of the holes. However, 
and as we discuss below, the determination of k F b,c is 
characterized by error bars in the order of ^ 1%, imply¬ 
ing an error of ~ 3% in the calculation of the volume of 
each pocket. Hence, the frequencies or cross sectional 
areas extracted at low fields are consistent with near 
perfect compensation, if one assumed that F a i and F a 2 
corresponded to hole pockets and Fp 1 and Fp 2 to elec¬ 
tron pockets. Instead, if one assumed for example, that 
two lowest frequencies, i.e. F a 1 and Fp± corresponded to 
the hole-pockets and the other two to electron pockets, 
which is the indexation suggested by Ref. [HI, one would 
obtain ^ 50 % more electrons than holes thus placing 
this system far from compensation. The other possible 
combination satisfying the indexation of Ref. [l5| would 
yield ^ 26% more carriers of one type with respect to 
the other. This would necessarily contradict the original 
claim in Refs. @,S of nearly perfect carrier compensa¬ 
tion. Hence, our indexation is far more consistent with 
the pronounced magnetoresistivity displayed by this com¬ 
pound. The FSs calculated by DFT, and subsequently 
adjusted to match the experimental results, lead to a 
much worse situation, with the density of electrons be¬ 
coming ~ 3 times larger than the density of holes. 


® 

<b) f - ; v 

<d) # V 

9 > € • 

FIG. 8. Evolution of the Fermi surface by shifting the Fermi 
level by 10 and 20 meV above and below the Fermi level, re¬ 
spectively. (a) The electron (shown in blue and transparent 
light blue color) and hole pockets (shown in red and in trans¬ 
parent pink color) for the actual position of the Fermi level. 
Here, the position of the Fermi level was adjusted to obtain 
a good agreement with the experimental Shubnikov de Haas 
frequencies for magnetic fields along the c-axis. (b) The posi¬ 
tion of the Fermi level is shifted upward by 10 meV. (c) The 
position of the Fermi level is shifted downward by 10 meV. 
(d) The position of the Fermi level is shifted upward by 20 
meV. (e) The position of the Fermi level is shifted downward 
by 20 meV. 

In Fig. [7] the geometry of the FS is compared to the 
geometry of the DFT calculated FSs. Figure 0 shows the 
position of the four fundamental peaks found if the FFT 
spectra of the SdH signal as a function of the angle rela¬ 
tive to the crystallographic c-axis. Here, 0° corresponds 
to the c-axis and ±90° corresponds to either the a- or the 
6-axis. Dashed lines depict the angular dependence of the 
frequencies resulting from our DFT calculations for both 
the hole- and the electron-pockets including spin-orbit 
coupling. These curves were obtained by a polynomial fit 
of their calculated values at discrete angles. Our calcula¬ 
tions yield a relatively good agreement with the experi¬ 
mentally obtained frequencies for fields along the c-axis, 
and a modest agreement with the angular dependence 
of the frequencies corresponding to the electron-pockets. 
Nevertheless, it shows a very poor agreement with the 
angular dependence associated with the hole-pockets. In 
fact the observed angular dependence imply very similar 
geometry for the electron- and the hole-pockets. Solid 
lines correspond to fits of the observed angular depen¬ 
dence assuming ellipsoidal FSs. As seen, and as already 
pointed out by Ref. 0 the experimental data can be 
well described in terms of ellipsoidal pockets. The cor¬ 
responding Fermi vectors k F b,c extracted from our fits 
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are given in Table II. Typical error bars in kp b,c are in 
the order of 1%. It is however quite remarkable that the 
DFT calculations fail to describe the electronic structure 
of this weakly correlated binary compound. 


F 

k a F (A- 1 ) 

k b F (A” 1 ) T 

kp (A 1 ) 

Vfs (A 3 ) 

F a 1 

0.04474 

0.06112 

0.12428 

0.0014237 

F a 2 

0.05668 

0.08578 

0.13886 

0.0020823 

Fpi 

0.05312 

0.07437 

0.12585 

0.0024678 

F/32 

0.05431 

0.08112 

0.13372 

0.0028278 


TABLE II. Fermi vectors and Fermi surface volumes resulting 
from the fits of the SdH frequencies as a function of the angle 
assuming ellipsoidal pockets. The error bars in the determi¬ 
nation of Kp b,c are in order of 1 %. 


In Fig. [ 8 ] we illustrate the effect on the geometry of 
the FSs of displacing the position of the Fermi-level by 
10 and 20 meV both above and below the Fermi surface. 
These displacements are intended to give us a rough idea 
of the role played by the Zeeman-effect on the geometry 
of the spin-orbit split FSs. The role of the Zeeman term 
is different for each one of the spin-orbit-coupling split 
pockets and different for electron and hole pockets and 
this effect is reversed when the magnetic field is applied in 
the opposite direction. Furthermore, the amount of the 
Zeeman splitting depends on k since the mixing of the 
atomic orbitals in forming the Bloch states changes for 
different values of k. Thus, we cannot infer the precise 
role of the Zeeman term by simply shifting the Fermi 
level as in Fig. [51 Nevertheless the results shown in Fig. [ 8 ] 
conveys the sensitivity of the geometry of Fermi surfaces 
sheets with respect to an applied field. Notice that the 
hole-pockets change significantly by displacing the Fermi 
level by a small amount in energy, i.e., of the order of just 
10 meV, which would correspond to a Zeeman splitting 
of the order of 100 T. The change is more dramatic when 
the Fermi level is shifted upwards. This exercise implies 
that the Zeeman-effect could induce pronounced changes 
in the geometry of the FSs altering the balance between 
opposite spin carriers. Notice also how the hole-pockets 
extracted from the DFT calculations, which present a 
boomerang like cross-section, nearly touch at its vertex 
when the Fermi level is displaced. Non-concentric Fermi 
surfaces like these would favor the magnetic breakdown 
between both orbits as experimentally seen. 

IV. SUMMARY 

In summary, by analyzing the angular dependence of 
the Shubnikov-de-Haas oscillations in single-crystals, we 
detect four Fermi surface cross-sectional areas in the 
semi-metal WTe 2 - By following density functional the¬ 
ory calculations, we conclude that they correspond to 
spin-orbit split electron- and hole-pockets or two sets of 


concentric pockets at either side of the T-point. Never¬ 
theless, the experimentally observed angular dependence 
for the cross-sectional areas is poorly described by the 
DFT calculations. Their geometry is well-described in 
terms of ellipsoidal pockets as suggested by Ref. 0. 
The discrepancy between the ab-initio calculations and 
the experimentally determined cross-sectional areas is far 
more pronounced for the hole-pockets. We observed four 
SdH frequencies (and not three as claimed by Ref. fl 2 h 
which are in excellent agreement with those reported in 
Refs. [H and El indicating that the geometry of the 
extracted Fermi surfaces does not dependent on sam¬ 
ple quality (for samples with resistivity ratios between 
10 2 < p(T = 300K )/p(T = 0.5K) < 2 x 10 3 ). We, as 
well as Ref El attribute the two middle frequencies to 
electron-pocket (s), with the other two corresponding to 
the hole one(s). 

In contrast to Ref El we attribute this geometry to 
the spin-orbit splitting of the hole and electron Fermi- 
surfaces. On the other hand Ref. [HI, based on their 
pressure dependent studies combined with band struc¬ 
ture calculations, attributes the middle two frequencies 
Fp i and F /32 to electron- and hole-orbits, respectively. 
Based on our angular-dependent study, this indexation 
would necessarily imply a sizeable difference in volumes 
between the electron- and the hole- Fermi surfaces, and 
therefore result in a pronounced lack of carrier compen¬ 
sation at zero-field. Here, we emphasize that our choice 
for the indexation of the Fermi surfaces in combination 
with the volumes of the extracted ellipsoidal Fermi sur¬ 
faces, do not discard near perfect compensation between 
electron- and hole-densities at zero- or low fields. 

By performing measurements up to higher fields we 
find experimentally that the effective masses in WTe 2 
are renormalized by the magnetic field, increasing con¬ 
siderably as the field increases. This is consistent with 
our fits of the oscillatory SdH signal as a function of the 
inverse field (to several Lifshitz-Kosevich components). 
These fits also indicate that the fundamental frequencies 
are evolving with field. Both observations point towards 
a magnetic field induced modification of the electronic 
structure at the Fermi level or a Fermi surface that is 
responsive to a magnetic field. Our DFT calculations as 
well as those in Ref. El indicate that the geometry of the 
Fermi surface in WTe 2 is very sensitive to the exact posi¬ 
tion of the Fermi level. Therefore, it is not surprising to 
find that the spin-up and spin-down Fermi surfaces can 
be slightly modified by the Zeeman-effect, which appar¬ 
ently is displacing the (spin-dependent) density of states 
at the Fermi-level towards van-Hove singularities. We 
believe that this effect could have a relevant influence on 
the observed colossal magnetoresistivity of WTe 2 . 

In a 3D Landau-Fermi liquid, the quasiparticle life¬ 
time decreases with increasing quasiparticle energy dif¬ 
ference from the Fermi energy [ 23 . If, upon application 
of a strong magnetic field the Fermi surface was altered, 
one should expect a strong influence on the quasiparti¬ 
cle scattering rates and lifetimes, contributing perhaps 
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in a non-trivial manner to the large magnetoresistance 
observed in WTe 2 - 

Finally, the lack of agreement between our experimen¬ 
tal observations and the band structure calculations has 
important implications. For example, the absence of a 
precise knowledge of the electronic structure at the Fermi 
level makes it difficult, or nearly impossible, to predict 
its topological properties 0- 
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